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Phenomenology of SUSY SU(5) with type I+III seesaw 

Carla Biggiqj and Lorenzo Calibblj 
Max-Planck-Institut fiir Physik (Werner-Heisenberg-Institut) , D-80805 Miinchen, Germany 

We consider a supersymmetric SU(5) model where two neutrino masses are obtained via a mixed 
type I+III seesaw mechanism induced by the component fields of a single SU(5) adjoint. We have 
analyzed the phenomenology of the model paying particular attention to flavour violating processes 
and dark matter relic density, assuming universal boundary conditions. We have found that, for 



a seesaw scale larger than 10 



GeV, BR(^ — >■ 67) is in the reach of the MEG experiment in 



sizable regions of the parameter space. On the other side, current bounds on it force BR(r — >■ fi'y) 
to be well below the reach of forthcoming experiments, rendering thus the model disprovable if a 
positive signal is found. The same bounds still allow for a sizable positive contribution to ek, while 
the CP violation in the Bs mixing turns out to be too small to account for the di-muon anomaly 
reported by the DO collaboration. Finally, the regions where the neutralino relic density is within the 
WMAP bounds can be strongly modified with respect to the constrained MSSM case. In particular, 
a peculiar coannihilation region, bounded from above, can be realized, which allows us to put an 
upper bound on the dark matter mass for certain set-ups of the parameters. 



I. INTRODUCTION 

Neutrino masses are an indication for the presence of 
new physics beyond the standard model (SM). The sim- 
plest extension consisting in adding to the SM fields three 
right-handed (RH) neutrinos and giving them a Dirac 
mass is not very satisfactory since it would require ex- 
tremely small Yukawa couplings, much smaller then the 
ones for charged particles. A nice way of explaining neu- 
trino masses as well as their smallness is through the so- 
called seesaw mechanism: new fields with masses much 
heavier than the electroweak (EW) scale, once integrated 
out, generate a dimension-five operator which, after EW 
symmetry breaking, gives a Majorana mass to neutrinos. 
When neutrino masses are generated at tree-level, three 
heavy mediators are possible, namely singlet fermions 
(corresponding to the so-called type-I seesaw [l|), triplet 
scalar (type H P) and triplet fermions (type III f3l|)"'^. 
Independently of the nature of the mediators, the neu- 
trino mass results to be suppressed by their heavy mass. 
Interestingly, 0(1) Yukawa couplings require the scale of 
new physics to be around the grand unification scale. It 
is then natural to study the seesaw mechanisms in the 
context of grand unified theories (GUT). Moreover, as 
it is well known, the presence of low scale supersymme- 
try (SUSY) triggers the unification of gauge couplings, 
so that usually SUSY GUTs are considered. 

In the literature plenty of models of SUSY GUTs in- 
cluding a seesaw mechanism for neutrino mass generation 
has been proposed. Actually, many different GUTs have 
been studied, involving the different seesaws, especially 
the types I and II. As for the type III (and type I-f-III), 
its embedding in a SUSY GUT was firstly proposed in 
Ref. Q and then discussed in Ref . 163 (in the context of a 



renormalizable model) and in Ref. [7|, and embedded in 
a fiavour model in Ref. Q- In all these cases the grand 
unified group considered is SU(5), which is the simplest 
group in which the SM gauge group can be embedded 
(for a wide discussion about SUSY SU(5) with the three 
different seesaw mechanisms, see 8]). In a minimal ver- 
sion of SU(5) [iO] neutrinos are massless, so that this 
GUT model has to be extended in order to account for 
neutrino masses. The model we consider here is some- 
how the simplest extension of a SUSY SU(5) accounting 
for neutrino masses, since with the simple addition of 
only one SU(5) representation two neutrino masses are 
generated via a mixed type I-t-III seesaw mechanism [7|. 
This is different from other models, where at least two 
representations were needed to account for two or three 
neutrino masses. Indeed, in the case of type I seesaw, 
the addition of at least two singlets is mandatory, while 
for the type II one must include both a triplet and its 
conjugate (or, in terms of SU(5) multiplets, a 15 and a 
15 [ll|)^. This model is the SUSY version of the model 
firstly proposed in Ref. [Ij, [3 ■ Also there two neutrino 
masses are obtained via a mixed type I-I-III seesaw, but 
being the model non-SUSY, the spectrum is completely 
different. Indeed, to guarantee unification, the triplet 
must be at the TeV scale, while the singlet mass is not 
specified. In this model lepton fiavour violation (LFV) is 
usually suppressed either by the small Yukawas or by a 
large mass. This is the standard situation in non-SUSY 
seesaw models, unless a low-scale inverse seesaw is real- 
ized 0, [Ij] . Furthermore, in presence of cancellations in 
the neutrino mass matrix, sizable rates of LFV processes 
are still possible (see for example, in the context of type 
I+III seesaw, [l5|b 

Here we study in detail the phenomenology of this 
"minimal" SUSY SU(5) with massive neutrinos, paying 
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^ See Ref. [j| for a review on the three mechanisms. 



^ Also in the type III models mentioned before more representa- 
tions are introduced, either matter [g| or Higgs [g]. 



particular attention to the effects on LFV processes. In 
order to isolate the effects purely induced by the RG run- 
ning of the SUSY parameters we assume universal soft 
masses at the GUT scale. We also discuss possible re- 
lated effects in the quark sector and how the region of 
the parameter space where the lightest SUSY particle 
(LSP) is a viable dark matter (DM) candidate are mod- 
ified. The rest of the paper is organized as follows: in 
Sect. In] we introduce the model, in Sect. UTT] we discuss 
the flavour violation, in Sect. IIVI we present numerical 
results and in Sect. |V] we conclude. The renormalization 
group equations (RGEs) of this model are gathered in 
the Appendix. 



II. THE MODEL 

We consider a SUSY SU(5) model where the matter 
content is enlarged with a 24 representation of SU(5), 
in order to get neutrino masses. Indeed the 24 can be 
decomposed under SU(3)^ x SU(2)^ x U{1)y asS + T + 
+ V + V, where 

S^{l,l)o r^(l,3)o 0^(8,1)0 U^(3,2)_5/6. 

S and T give mass to two neutrinos via a mixed type 
I+III seesaw mechanism. The relevant superpotential 
terms are then 

Wsu{5) - ^(Yio)y 10aO,5H + \/2(r5)», 10,5,5h + 

2/^4 5,24 5h + ^24 24 + ••• (1) 

Here 5 = (d^L), 10 = (u^Q,e'=) and 5h = {Ht,Hi), 
as usual. Below the SU(5) breaking scale, the superpo- 
tential reads: 



the resulting misalignment is too small, giving rise to a 
second neutrino mass much smaller that the solar mass. 
It is then clear that the GUT relation on the Yukawa 
couplings must be somehow altered, in order to get two 
massive neutrinos. Since non-renormalizable operators 
should anyway be added to the SUSY SU(5) lagrangian 
in order to break the unwanted relation l^e = Yj which 
is in disagreement with the experimental value of fermion 
masses [16['^, we can also add non-renormalizable opera- 
tors involving the 24 jl3 |: 



W'nr 3 \ [y\ (5, 24)(24h 5h) + vl (5. 24h)(24 5h)- 



y^5,tr(24 24H)5H] 



(4) 



where A represents the cut-off of the theory (e.g. 
the Planck scale M-pi) and 24h is the Higgs in 
the adjoint representation which breaks SU(5), with 

(24h) = diag(l/3, 1/3, 1/3, -1/2, -1/2)^/3/5 WG. The 
three above operators could be for instance generated by 
integrating-out, respectively: a Higgs superfield in the 
45-1-45 representation, a matter 45-f 45, a SU(5) singlet, 
all having a mass of the order of the scale A. 

After SU(5) breaking, the Yukawa couplings of the dif- 
ferent components of the 24 split and the matching at the 
GUT scale reads: 
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W = {Yu)^Ju1QjK + {Yd)^jd1QJhd + {Y^),Je1LJhd + 
2/y LiThu + Vs LiShu + yy dlVK + 

M±SS + ^TT + ^00 + My VV. (2) 

AAA 

Notice that, from Eq. ([T]), it follows that the S'C/(3) octet 
O and the field V do not have Yukawa interactions with 
fields which remain lighter than the GUT scale. 

From Eq. ([2]), it is easy to see that the singlet S and 
the neutral component of the triplet T generate neutrino 
masses via a seesaw mechanism: 



!i (ysyi.yryr^ 

2 y Ms Mt j 



(3) 



Notice that, since S and T belong to the same SU(5) 
multiplet, at the GUT scale yr — ^10/3 ys — 2/24 /"\/2. 
The previous rrii, is then a rank-1 matrix and only one 
neutrino mass is generated. In principle the Yukawas 
could be misaligned in the running from the GUT scale 
to the seesaw scale simply via RGE effects, but in practice 



In practice the coupling 2/3 determines the misalignment 
between ys and yV that permits to generate two non-zero 
neutrino masses. As we will see later, LFV processes 
suggest for the seesaw scale an upper bound around 
2Q13-M4 QqY, corresponding to neutrino Yukawa cou- 
plings smaller than 10~^^^^'. From Eq. ^ we then de- 
duce that 1/3 should be 0(1) and the cutoff scale A cannot 
be too large, if we want to saturate these bounds. From 
now on we will consider ys and yr as independent pa- 
rameters. 

In the same way as the splitting among the Yuwakas 
is generated via SU(5) breaking effect, also a splitting in 
the masses arise. The terms contributing to the masses 
of the 24 fields are 



1 /v^ 

-M24tr (242) + A24tr (24^ 24h) + ^ ( ^ 



(8) 



^ An alternative solution is adding new Higgs representation [17|| . 



giving rise to the following masses: 
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(10) 

(11) 
(12) 



where the actual expressions of the A^^^ in terms of the 
possible NR operators couplings can be found in |12| . 

If we do not allow fine cancellations, we can conclude 
that the components of the 24 have masses of the same 
order of magnitude, so that only small threshold effects 
are introduced and the successful gauge coupling unifica- 
tion of the MSSM is maintained. Notice that, since LFV 
processes disfavor a seesaw scale larger than lO^"* GeV, 
the coupling A24 should be 10"^^"^^-' or smaller while the 
non-renormalizable couplings Xf^ can be larger, depend- 
ing on the cutoff scale A. 

We now focus on neutrino masses. Eq. ^ can be in- 
verted and the Yukawa couplings can be expressed in 
terms of low energy parameters. Since in this case we 
have only two massive neutrinos, Casas-Ibarra parame- 
terization is simplified, since the R matrix depends now 
only on one complex parameter z |18l| . The Yukawa cou- 
plings can then be expressed as |13j : 



vy¥ 

and 




vys 




H3^ 



m\ sin z — Ui2 a 



'rrig cos z) NH 
/rn^cosz) IH 

(14) 
There is another solution with the opposite sign for the 
second terms in Eqs. ((T^ [H)) . In the above equations, 
Uij are elements of the PMNS matrix^. For the neutrino 
mass eigenvalues we have in the normal hierarchy (NH) 



m^ 



where we take the neutrino mass parameters [21| as mea- 
sured in the solar and atmospheric oscillation experi- 
ments 



Ato| « 7.59 X 10"^ eV^ 



Ar 



(17) 
2.46 (2.36) X 10"^ eV^ NH (IH). (18) 



Ami mg = V Am^ + Am|, (15) 



From Eqs. (IT^ [HI) we see that the Yukawa couplings 
grow with the square root of the mass of the heavy 
fermions, as a trivial consequence of the seesaw formula 
Eq. ([3]), and increase exponentially when Im(z) grows. If 
we want to avoid unnatural cancellations in the neutrino 
sector (e.g. between the two terms of Eq. (jH])), Im(z) 
should be < 0(1). As we will show later, the present 
bound on /i — > 67 can actually constrain it to smaller 
values. 

Notice that in what respect neutrino masses, this 
model is not different from a model with two right- 
handed neutrinos (2RHN) |18l - l20| , where for instance the 
same parameterization of Eqs. p^ [T4l) holds. However, 
this model, besides the fact of being better motivated 
from a GUT perspective, has got some features which 
distinguish it from a generic 2RHN model: 

• Up to fine-tuning, the parameter space is more re- 
stricted, since it is natural to assume Ms and Mt to 
be of the same order of magnitude (see Eqs. (P ITU)) ). 
Moreover, barring cancellations, y's and y^ will be 
also of the same order of magnitude (Eqs. (0 16])). 

• The presence of the SU(5) partners of S and T in- 
duces fiavour violating effects in the hadronic sector 
as well, similarly to what happens in the leptonic 
sector. Again, the relevant couplings j/y, even if 
in general independent, are expected to be of the 
same order of y's and y^. 

• The presence of a full 24 at an intermediate scale 
between the GUT and the EW scales does not spoil 
gauge coupling unification if Mt — Mq — My — 
All , as in our case, but affects the gauge couplings 
running above M/ . This can have an impact on the 
SUSY spectrum and, in particular, on the regions 
of the parameter space which provide a relic density 
for the LSP within the WMAP bounds, as we are 
going to discuss in section ITVl 



while in the inverted hierarchy (IH) case neutrino masses 
are given by: 



TOi = \/Am^ - Ato| m^ = J ^m\ TO3 = 0, (16) 



* We remind the reader that here, hkc in the 2 RH neutrinos case, 
the PMNS matrix only has two phases: a Dirac phase 5 and a 
Majorana phase <&. 



III. FLAVOUR VIOLATING PROCESSES 

The presence of the fields in the 24 modifies the renor- 
malization group running of the parameters of the model, 
both the superpotential couplings and the SUSY break- 
ing terms. For instance, the renormalization group equa- 
tions (RGEs) for the scalar masses are now given by: 



1G 2 ^ 2 qMSSM,! 

16- J^rn, = /?,„. 



B^^^\ 



(19) 



where /3 2 'is the usual MSSM 1-loop /3-function and 



r 



is the new 1-loop contribution given by the new 



fields in the 24, with clearly /3 2' 7^ only above the 24 

energy scale. 

In particular, the couplings of the seesaw fields, S and 
T, with the lepton doublet will affect the running of the 
left-handed (LH) slepton masses, generating off-diagonal 
flavour violating entries, in perfect analogy with what 
happens in the context of supersymmetric seesaw of type 
I [l^, 122] and type II [ll| . In addition, the presence of 
the SU(5) partner, V , of the seesaw fields will induce an 
analogous effect for the RH down squarks. 

The complete RGEs of the model are given in the Ap- 
pendix. Let us display here the /3-functions of, respec- 
tively, the LH slepton and RH down-squark soft masses, 
which are the relevant ones for outlining the effects men- 
tioned above: 



L ' IJ 



2 Vn Vt, (ml +m%) + 2 A*^^, A. 
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ykivs^Vj + imly*shysj + 
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2 yh yh (™k + ™|) + 2 A*s, Al^) , (20) 

^It^j = ^yvt (yv m%)j + {m\ yy), y^^ + 

d^ / ij ^ 

2 yh yvj ("^k + "^1) + 2 A*y, A^,) .(21) 



Off-diagonal flavour violating entries in the LH slepton 
and RH down-squark mass matrices are then generated 
by RG running from Mqut down to the seesaw fields 
mass scales, even starting with universal boundary con- 
ditions at Mgut, m| -^ ^ mil. From Eqs. (HOI HH), 



we can estimate the flavour violating mass-insertions 



mZiTi^^, which parameterize the amount 



of flavour violation induced by the running. At leading- 
log, they read: 
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(22) 
(23) 



where to^ , m|^ are average slepton and squark squared 
masses at low energy. Eqs. ([^ [^ provide a good esti- 
mate of the FV mass-insertions, unless ttiq is too small. 
In the case of mo ~ 0, which is indeed possible in the 
model as we will discuss in the next section, the sfermion 
masses are generated by the running, but Eqs. (1^ [^5)1 
are clearly not valid anymore, since sfermion masses are 



vanishing at AfcuT and possible off-diagonal entries in 
the mass matrices can be only generated at orders higher 
than the leading-log. 

Keeping that in mind, we can still make use of Eq. (P^ 
to get an idea of the expected amount of LEV. Eor in- 
stance, in the case Ms — Mr — 10^'^ GeV, we see 
from the seesaw formula, Eq. ^ , that typically y^ j, ~ 
0(10"^) and therefore assuming m~ 
Eq. ([^ gives roughly: 



TTiQ and ^0 — 0, 



(<5£l),j ~ 0(10-*) for Ms^T ^ 10^^ GeV 



(24) 



value which can give sizable effects in the /i-e transitions 
only and can already exclude the SUSY parameter space 
in the light sleptons regime^. 

We can also estimate the typical ratio of the BRs of 
different LEV processes. Given that the main source of 
LEV is represented by the {Sf^i^)ij, we have: 



BR(£, ^ £,j) 



« \{5li.h? 



(25) 



hence the ratio of BRs in the t-^ and ^-e channels can 
be estimated to be: 



R 



BR(t -^ n-f) 
BR(^ -^ 67) 
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(26) 



Using Eqs. (II31II11), one can check that 4 < i? < 80 in the 
normal hierarchy case for a real parameter z, if Ues — 
and the Majorana phase $ vanishes as well. As the value 
of Ue3 increases, one can verify that R diminishes and 
it becomes 0(1) for Ues — 0.2. Interestingly, as soon 
as Im(z) is switched on, R rapidly drops to 0(1) values 
independently of the value of Ues ■ As for the role of the 
phases, they also generically tend to reduce R, even if for 
small non-zero values of C/gs the Dirac phase somehow 
compensates the Ues effect, preventing the reduction of 
the ratio. Moreover, the presence of the phases increases 
((5£l)i2, from which a bound on Im(z) can be derived (see 
later). In the inverted hierarchy case with a real z, R 
can even diverge, since for certain values of z (i5£l)i2 can 
vanish. However, such cases correspond to set-up of the 
Yukawas (e.g. y^j, ^ J/ts) which cannot be considered 
natural in the light of Eqs. ([5115]). Moreover, for Im(z) ^ 
0, such divergences disappear and R tends to 0(1) values 
like in the case of normal hierarchy. 

Let us briefly make here a comparison with other see- 
saws, still implemented in a SU(5) context. As already 
discussed in the previous section, in what respect neu- 
trino masses our model is not different from a model with 
a type I seesaw with only two RH neutrinos. This state- 
ment holds also for LEV, with the only difference given 
by the fact that in the 2RHN model the heavy neutrino 



Cfr. the bounds provided in [23 



masses can be hierarchical, while our model, barring can- 
cellations, predicts Mx — Ms- As a consequence higher 
values for R can be obtained [20|. In the type I seesaw 
with three RH neutrinos, due to the larger number of pa- 
rameters, even more freedom is allowed. On the contrary 
in the type II seesaw there is a direct relation between 
high-energy and low-energy neutrino parameters, so that 
the ratios of the branching ratios of LFV processes can 
be expressed in terms of neutrino masses and mixing an- 
gles. When this is embedded into a SU(5) GUT by adding 
a 15-1-15 representation [ll|, R varies between 400 and 
0(1) with increasing 6*13 |24| . Notice that in that model 
the seesaw fields induce flavour violation in the hadronic 
sector too, as in the model we are studying in this paper. 
Let us now discuss the induced flavour violation in the 
hadronic sector. From Eq. (P5|) . we see that off-diagonal 
entries in the rn^ mass matrix are induced by the cou- 
pling yv of the down-quark SU(2) singlets with the 24 
field V. Eqs. dSJT]) tell us that yv cannot be unequivo- 
cally determined in terms of ys and j/r and, therefore, in 
terms of neutrino parameters. However, Eqs. (IS][7]) also 
show that yv can be naturally expected to be of the same 
order of magnitude of the seesaw Yukawa couplings, as it 
clearly follows from the SU(5) embedding of the model. 
In particular if |27i,2| < \y24\ m Eqs. (|5T7|) . then yv ~ yr- 
In our numerical analysis, we are going to make use of 
this last assumption^. Anyway, with a free choice of 2/1,2 
and barring cancellations, we would still get: 



i^RRjij 
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(27) 



Thus, comparing this expression with Eq. (|24p . we find 
that the typical order of magnitude of the hadronic mass 
insertion is: 



(^i 



rr)u 



0(10"^ H- 10"-*) for Mv ^ 10'^ GeV, (28) 



with ((5pf{)ij becoming maximal for toq 3> Afi/2, when 
mj ~ nil — rrin. Moreover, one has to take into ac- 

d^ L ^ ' 

count that, like in the MSSM, the RGE for m^ gener- 

ates also small {Sf^i^)ij proportional to CKM elements, 
('^ll)«j ^ ^ti^tj- Taking into account this further effect, 
the most stringent bounds, which come from the Kaon 
system, are |(5rr)i2| — O(10~^) [11]. Comparing this 
value with Eq. ((28)) . we see that the model does not typi- 
cally predict large deviations from the constrained MSSM 
(CMSSM) predictions in the hadronic sector and there- 
fore it is safe from hadronic FCNC constraints. However, 
Eq. (j28p provides a quite rough estimate and it is worth 
to study in more detail some hadronic observables, for 



Another possible approach to improve the predictivity of the 
model in the hadronic sector is considering the renormalizable 
version of the model discussed in [g| , where yv can be written as 
a combination of yx and ys- Apart from this point, that model 
gives the same phenomenology discussed here. 



which experiments have recently showed possible ten- 
sions with the SM (and the CMSSM) predictions. In 
the next section, in particular, we are going to comment 
about the impact of the new flavour mixing sources of 
the model on the Kaon CP-violating parameter ex and 
on the time-dependent CP asymmetry, S^^, in the decay 



IV. RESULTS 

As mentioned above, in order to outline the effects in- 
duced by the RG running between the GUT scale and 
the mass scale of the 24 fields, we consider universal 
boundary conditions, namely: a universal scalar mass 
Too, a common gaugino mass Mi 12 and trilinear terms 

In order to compute the SUSY spectrum, we numer- 
ically solve the full 1-loop RGEs of the model (see the 
Appendix) down to the seesaw scale, Als = Mt = Mj, 
at which the 24 fields decouple. Then, we run the MSSM 
RGEs down to the SUSY scale msusY = ^/^Wj^- Fo^' 
each point of the parameter space, we impose the fol- 
lowing requirements: (i) successful EWSB and absence 
of tachyonic particles; (ii) limits on SUSY masses from 
direct searches; (iii) neutral LSP. Then, we compute 
the leptonic processes by means of a full calculation in 
the mass eigenstate basis [2^, the hadronic processes by 
means of the mass-insertion approximation formulae in 
[27| , the LSP relic density using DarkSUSY [28] and the 
BR{B -^ Xsj) using SusyBSG ^. We require that the 
resulting BR{B — )■ Xgj) do not deviate from the experi- 
mental value [30[ in more than 3cr. 

Let us first try to extract information about the seesaw 
scale and the other seesaw parameters, focusing on the 
stringent bounds LFV can impose on them. In order 
to do that, we can start varying all the parameters in 
large ranges, but we clearly need a criterion for defining 
the SUSY spectrum we want to concentrate on (all effects 
would be negligibly small, if we considered slepton masses 
of several TeV). Therefore, we will mostly concentrate on 
parameter regions giving sizable SUSY contributions to 
the anomalous magnetic moment of the muon, (g — 2)^, 
and, later on, also on regions which provide a dark matter 
relic density within the WMAP constraints. 

For simplicity, in the numerical analysis we have ne- 
glected possible 0(1) mass-splittings among the fields 
in the 24, inducing threshold corrections to the gauge 
coupling running, which would modify the MSSM gauge 
coupling unification. In particular, the 1-loop prediction 
for the value of the strong coupling at AIz would become: 
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where a'^{AIz) is the 1-loop MSSM prediction. According 
to the above formula, the consistency with the measured 
value for 013 (M^) could be slightly worsened or improved. 
We notice that such modification of the running of the 
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FIG. 1: BR(/i -^ 67) vs. the seesaw scale Mi (in GeV) for 
tan /3 = 10 (top) and tan /3 = 40 (bottom) and a wide scan of 
the parameters (see the text for details). Blue points provide 
(5-2)fSY/2>io-^ 



gauge couplings would have anyway a negligible impact 
on the running of the other parameters. Moreover, the 
possible thresholds would affect the running of the soft 
masses, entering in the expressions of the flavour violat- 
ing parameters, Eqs. (|22l [23|) . only logarithmically, so 
that they would have a small effect on the observables 
we are going to study. 



A. Lepton Flavour Violation 

In Fig. [1] we plot BR(/Lt — >■ 67) as a function of the 
seesaw scale Ms = Mt = M/, in the case of normal 
neutrino hierarchy^, for the following choice of the pa- 
rameters: tan/3 = 10 (top panel) and tan/3 = 40 (bot- 
tom panel), < too < 1 TeV, < M1/2 < 1.5 TeV, 
1^0 1 ^ Stoq. The neutrino parameters were also varied 
in the following ranges: 10^° GeV < Mj < 10^^ GeV, 
10-3 < \Ue3\ < 0.2, < z <2tt. We took the parameter 



''' We checked that inverted hierarchical neutrinos do not provide 
significantly difi'erent predictions with respect to the normal hi- 
erarchy case. We will thus concentrate on normal hierarchy from 




FIG. 2: Contour plot for the LFV parameter |((5£l)i2| in the 
Re(2)-Im(2) plane, for mo = M1/2, Aq = 0, Mi = 10^^ GeV 
and Ue3 = 0, $ = 0. 



z real, since the only effect of its imaginary part is to raise 
the seesaw Yukawas and so the /i — >■ 67 rate, as we com- 
mented in Sec. HIl However, BR(/i — )■ ej) itself provides 
very stringent bounds on Im(z), as we will comment be- 
low. We have also checked that all couplings remain per- 
turbative up to the GUT scale. The blue (black) points 
give a^usY = (g - 2)^usY/2 > iQ-s, so lowering the 
tension between theoretical prediction and experiments 
below the 2a level. 

For tan/3 — 10 (upper panel of Fig. [T]), we see, besides 
the dependence BR(/i — > 67) ^ Mf, that the current 
experimental limit BR(/i -^ 67) < 1.2 x 10"-^^ j3ll|, con- 
strains the seesaw scale to be Mj < 10^^ -^ 10^** GeV 
for the points favored by {g — 2)^, even if there are 
few points, for which the parameters conspire in low- 
ering BR(/i — > 67), that can evade such bound. Even 
if BR(/i — >■ 67) is enhanced by increasing tan/3, we find 
the above limit on the seesaw scale also for tan/3 = 40 



(lower panel of Fig. [ij , since a^^^^ is increased by tan /3 
as well. In both cases, the MEG experiment [S^l, whose 
expected sensitivity is BR(/i — > 67) ~ 10-^^, will be able 
to test soon the region of the parameter space favored by 
{g - 2)f, down to M/ ~ lO^^ ^ lO^^ GeV. 

Let us now show how the present bound on BR(/i — > 
e"f) can severely constrain the parameter z. We have al- 
ready argued above that Im(z) cannot be too large, with- 
out having unnatural cancellations in the neutrino mass 
matrix m,^ between the singlet and the triplet terms. 
Besides that, LFV bounds can directly constrain Im(z), 
since the seesaw Yukawas simply grow by increasing it. 
For convenience, let us express the BR(^ -^ ej) con- 
straints in terms of bounds on the mass insertion (^£l)i2- 
From the same scan of the parameters presented above, 
we find that satisfying the present limit on /i — > 67 re- 
quires: 



|(^£l)i2|<(5-6)x10- 



(30) 
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FIG. 3: BR(^ -> 67) vs. BR(r -^ fi-y), for tan/3 = 10 and the 
same scan of the parameters of Fig. [l] 



for the points lying in the blue (black) region. In Fig. [21 
we show contours for |(<5£l)i2| in the Re(z)-Im(z) plane, 
for Too = Ml/2, ^0 = 0, Mj = lO^^ Q^y ^nd Ues = 0, 
$ = 0. We see that, indeed, |(^ll)i2| grows very fast 
with Iin{z). As a consequence, the bound of Eq. ([50)1 
constrains Im(z) to values ^ 0.3 for the seesaw scale at 
10" GeV. 

Let us finally consider LFV in the fi-T sector as well. In 
Fig.m we plot BR(^ -> 67) vs. BR(t -> ^7), for tan/3 = 
10 and the same variation of the other parameters of 
Fig. [TJ We see that in this model, the strong bound 
on flavour transition in the /i-e sector already challenges 
future T — >■ /X7 experiments quite strongly. In fact, the 
bulk of the points, for which BR(/i —> e-f) is less than 
the present bound, gives BR(r — > fj,j) ^ 10~^, which 
is indeed below the expected sensitivity of the proposed 
Super Flavour Factory (Sj). This is consistent with the 
estimate for R = BR(r -J> ^7)/BR(/i -)> ej) < 100, 
we provided in the previous section. However, we see 
that there are some points for which parameters conspire 
to raise the value of BR(t — )■ fij) at the level of 10 
i.e. in the reach of the SuperB factory at KEK 



M\ 



Nevertheless, a positive signal for BR(r 
anyway disfavor the scenario under study. 



fij) would 



B. Hadronic observables 

In this section, we are going to discuss the effects of 
the new source of flavour mixing in the down squark sec- 
tor, Eq. (1231), induced by the running between the GUT 
and the seesaw scales. In particular, it is interesting to 
check if this is able to account for a large phase in the Bg 
mixing, as suggested by the Tevatron experiments CDF 
|35| and DO [3^, |33|- Moreover, a positive new physics 
contribution to ex (around the 24% of the SM contribu- 
tion) [3^ is one of the possible ways for accommodating 
a recently reported tension among different observables 
used to fit the unitarity triangle (see also [2^ ,39] ) . 

As pointed out in section Hill hadronic flavour mixing 
cannot be directly related to the leptonic one. Never- 
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FIG. 4: Top: BR(^ -> 67) vs. ex/e^^ for 10 < tan/3 < 
40 and the same scan of parameters as in Fig. [H Bottom: 
BR(^ -^ e-y) vs. S^^. 



theless, the off-diagonal entries of m'^^ are naturally of 
the same order of magnitude as the leptonic ones. For 
definitiveness, we take yv = yr as input for the RGEs at 
the GUT scale, then we vary the phases of the resulting 
low-energy mass insertions {5^^)ij between and 27r. 

Let us start to consider the possible effect of the gen- 
erated ((5Rf{)i2 on the SUSY contribution to e^- In the 
top panel of Fig. m we plot BR(/i — >• 67) vs. ck/ea^ 
for 10 < tan/3 < 40 and the same variation of the pa- 
rameters as in Fig. [H As for the previous plots, the blue 
(black) points provide a sizable SUSY contribution to 
{g — 2)^, while the green (light grey) points give a neu- 
tralino relic density not larger than the cold dark mat- 
ter relic density measured by WMAP (see the next sec- 
tion for details). We can see that the present bound on 
BR(/i — >• 67) still allows for a sizable (up to 20-;-30 % of 
^k) positive contribution to €k- Furthermore we have 
BR(^ -^ 67) > 10"" (so within the sensitivity of MEG) 
for most of the points which provide such a solution to 
the ex tension, which would be therefore strongly disfa- 
vored by a negative result of MEG. Let us notice, how- 
ever, that the parameter space points favored by WMAP 
cannot provide the desired increase of ck- The reason is 
that these points are mainly concentrated in the coanni- 
hilation region where M1/2 > thq or even M1/2 ^ "T-Oj 
as we are going to discuss in the next section, so that 
the flavour violating {5^-^ij result suppressed by large 
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squark masses. 

In the bottom panel of Fig. SI we plot 5*^0 for the 
same scan of the parameters. As we can see, the pre- 
dicted value do not deviate too much from the small SM 
prediction S*^^ ~ 0.036. The reason is that, even if the 

phase of ((5ppj)23 can be large, |(^rr)23| is numerically 
too small (cfr. for instance the estimate in Eq. (1^51) ) to 
provide a sizable CP violation in the Bg mixing and thus 
accounting for the di-muon anomaly reported in [37|. If 
such new physics effects in Bg mixing will be confirmed, 
the minimal version of the model we are discussing here 
should be extended to include further sources of flavour 
violation in the hadronic sector. 
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C. Neutralino relic density 

The presence of intermediate scale fields, which are 
charged under the SM gauge group, has a possible im- 
pact on the supersymmetric spectrum and, thus, on the 
parameter space regions, for which the relic density of 
the LSP (in our case a bino-like lightest neutralino as 
in the CMSSM) resuhs to be within the WMAP bounds 
|40| . In this section, we are going to focus on the so- 
called f coannihilation region [41|, since focus point 42 1 
and A-funnel [43 are not expected to be qualitatively 
different with respect to the CMSSM (even if they can 
be quantitatively modified, even significantly, within this 
model) . 

The effect we are going to discuss can be again traced 
back to the modification of the RG running of the param- 
eters. In this case, however, this is not due to the new 
Yukawa interactions (since fiavour bounds do not allow 
the couplings to be too large), but it is an effect of the 
modification of the running of the gauge couplings (and 
the gaugino masses) above the scale of the 24. In fact, 
even if the fields in the 24 do not spoil (at least at 1-loop) 
the successful gauge coupling unification of the MSSM, 
the running gets "stronger": above Mt — My — Mq 
the 1-loop /3-function coefficients gets indeed modified as 
follows: 



-MSSM 



h-24 



(33/5, 1,-3) + (5, 5, 5), (31) 



and the running of the gauge couplings is considerably 
deflected. As a consequence, even if the couplings unify 
at the usual MSSM GUT scale, Mgut ~ 10^^ GeV, the 
value of the unified coupling au gets larger than in the 
CMSSM. 

Clearly, an analogous effect happens to the gaugino 
masses, so that they reach values at Mj, which can be 
considerably smaller than the unified value Mi/2. This 
could be thought as a simple rescaling of M1/2 (since 
clearly the low-energy gaugino masses will be the same 
as in the CMSSM with a lower value of M1/2), if it did 
not affect the scalar masses as well. In fact, with the 
same values of the gaugino masses at low energy, the 
scalar mass will feel a stronger gauge contribution to the 



FIG. 5: Coannihilation region and BR(/^ — > ej) contours in 
the (mo, M1/2) plane, for tan/3 = 10, Mi = 10^^ GeV. 



running, through the gauge terms in the RGEs, ~ aiMf , 
which are larger than in the CMSSM between the GUT 
scale and Mj. The consequence in the low-energy SUSY 
spectrum is that the scalar masses will result relatively 
larger, with respect to the gaugino masses, than in the 
MSSM. 

Qualitatively, the above described effect is clearly com- 
mon to all models that have fields charged under the SM 
gauge group at some intermediate scale, and it was, for 
instance, observed in the context of an 5*0(10) type-II 
seesaw model in 44^ and in a multi-scale fiavour model 
in 0. 

Coming back to DM, having relatively heavier scalars 
could destabilize the ordinary regions of the parame- 
ter space that provide a neutralino relic density within 
the WMAP bounds and for which quite precise relations 
among parameters are usually required. An example is 
the f coannihilation region, where the correct relic den- 
sity is achieved thanks to an efficient f-x? coannihilation, 
1^0 . As we are going to see, such 



which requires nif-^ 

region is strongly modified in our case, as an effect of 
the relatively heavier fi resulting from the strong gauge 
running below the GUT scale. 

What can happen to the coannihilation region is de- 
picted in Fig. [5l where we show the mQ-Mii2 plane for 
Mi = 10" GeV, tan/3 = 10, Aq = 0. The neutrino pa- 
rameters (not relevant for DM) were taken to be C/e3 = 0, 
z — 1. The region marked with red (grey) points gives 
^DMh"^ < 0.13. We can see that the CMSSM region 
where rrif-^ < m^o, along which usually the coannihi- 
lation strip runs, has disappeared as a consequence of 
the effect described above*. Coannihilation is still pos- 



* This opens up the possibihty of having mo = 0, i.e. vanishing 
scalar masses at high-energy (then generated through the run- 
ning driven by the gaugino masses), such as in [44L l45t . This 
possibility has been recently addressed in 1491 . 
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the value of ajj, and vice- versa would become stronger 
for lower values of Mj. A larger value of tan 13 would con- 
tribute as usual to decrease the fi mass (by increasing 
the negative contributions oc y^ in the running and also 
the L-R mixing term in the f mass matrix). This can 
be seen in Fig. [BJ where the tan /3 = 40 case is shown. 
The parameter space is now qualitatively similar to the 
CMSSM case: the region where fi is the LSP has reap- 
peared and the coannihilation region is a strip along it. 
Notice the presence for low values of mo and Mi/ 2 of a 
sizable "bulk" region (which is smaller but still present 
also for tan/3 — 10). This region is, however, already ex- 
cluded by several constraints, including the experimental 
limit on BR(/x — >■ ej). 



FIG. 6: Coannihilation region and BR(/^ — >■ ej) contours in 
the (mo, M1/2) plane, for tan/3 = 40, Mi = 10^^ GeV. 



sible, since very low values of mo still gives m^-^ ~ rn^o , 
but, interestingly, such region is bounded from above: 
this means that this particular set-up of the parameters 
predicts an upper bound on the DM mass, in this case 
TO^o ^ 240 GeV, as we can see from the figure taking 
into account that, for Mj = 10^"^ GeV, the bino mass is 
approximately Mi ~ 0.31 Mi/2- A similar effect, provid- 
ing an upper bound on m j^o , was found in [43] , again as 
a consequence of the modification of the gauge contribu- 
tion to the running of the scalar masses^ (in that case 
an additional SU(5) running of the parameters above the 
GUT scale was taken into account). 

In Fig. [51 we have also plotted contours for BR(/x — > 
67), the LEP limit on the Higgs mass (taking into account 
3 GeV of theoretical error), as well as the region which 
provides {g — 2)^^^^/2 > 10~^ (below the magenta dot- 
dashed line). We can see that the DM region is already 
partially excluded by the present limits on /i — >■ 67 and 
on the Higgs mass. The rest of the coannihilation region, 
which is, at least in part, consistent with a sizable {g — 
2^SUSY^ will be fully tested very soon by MEG, since 
it gives BR(^ — > ej) > 10^^^. The prediction for the 
fi ^ e"f rate clearly depends on the parameter z, which 
we have here taken z = 1. Nevertheless, we checked that, 
varying the value of z, still BR(/i — ?> 67) is predicted in 
the reach of the MEG experiment (i.e. > 10~^^) in the 
coannihilation region, apart from few points where the 
combination of parameters happens to suppress the rate. 

Finally, it is important to stress how the effect de- 
scribed above and its possible impact on the coannihila- 
tion region are sensitive to variations of the parameters, 
especially Mj and tan/3. The effect would be clearly de- 
creased, and would eventually disappear, by increasing 
Mi, i.e. decreasing the length of the running and thus 



See also \m]5> 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have considered a SUSY SU(5) model 
where neutrino masses are obtained via a mixed type 
I-I-III seesaw mechanism and we have studied its phe- 
nomenology assuming universal soft masses at the GUT 
scale. The main characteristic of the model is the pres- 
ence of one massless neutrino. Then the high-energy see- 
saw parameters are less than in the three massive neutri- 
nos case and therefore a higher degree of predictability 
is present. Moreover, the model represents a very eco- 
nomical way of accounting for neutrino masses in a GUT 
context, since the addition of just one chiral superfield in 
the SU(5) adjoint representation is considered. 

Besides discussing the model and in particular the 
mechanism through which we obtain two neutrino 
masses, we have analyzed the following features: 

• /x — >■ 67 and other LFV processes; 

• possible contributions to hadronic observables; 

• neutralino relic density. 

We have shown that we can have sizable contribu- 
tion to the /i — >■ 67 rate, such that the current exper- 
imental limit constrains the seesaw scale to be Mi < 
10" -=- 1014 GeV, while MEG wiU be able to test the 
model down to scales of Mj < lO^^ ^ xq" GeV. We 
have also shown that the bounds on BR(/i -^ 67) put 
strong constraints on BR(t — > fj,'-/) , making very unlikely 
to observe it in future experiments. Otherwise, a positive 
signal for the r — ^ fi^ decay would disfavor this model. 

From the bound on BR(/i — > 67), we have been able to 
put an upper bound on Im(z) ^ 0.3, for the seesaw scale 
Mi — lO^'^ GeV. Of course this bound is Af /-dependent, 
since a reduction in the scale would imply a decrease 
of the size of Yukawas and then larger values for Im(2:) 
would be allowed. However, as discussed in Sect. [Hi val- 
ues of Im(2:) larger than 1 are unnatural since cancella- 
tions in the neutrino sector would be needed. 

The contribution in the hadronic sector is given by the 
coupling of the down-quark singlets with the new fields 
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V. Even if this cannot be directly related to the neu- 
trino parameters, an order of magnitude estimate can 
be performed. We have shown that in this model the 
present bound on BR(/j, — >■ ej) still allows for a siz- 
able (up to 20-^30 % of e^) positive contribution to ek, 
which would help in accommodating a recently reported 
tension among different observables used to fit the uni- 
tarity triangle. On the other side, CP violation in the Bg 
mixing turns out to be too small to be able to account for 
the di-muon anomaly reported by the DO collaboration. 

As for the neutralino relic density, we have focussed 
on the so-called t coannihilation region. We have shown 
that the CMSSM region where m^^ < m^o, along which 
usually the coannihilation strip runs, can disappear, so 
that the coannihilation region gets distorted. Interest- 
ingly, such region is bounded from above, which means 
that an upper bound on the DM mass can be derived: 
for the particular set-up of the parameters considered 
here we got m^o ^ 240 GeV. Moreover, the possibility of 
having ttiq = (with efficient coannihilation as well) as 
high-energy boundary condition is now open. 

In this paper we have not addressed other issues such 
as proton decay and leptogenesis. As for proton decay, 
this model does not improve the situation with respect to 
the standard case, so one has to rely on standard mech- 
anisms to suppress the proton decay rate. For instance, 
in the context of the missing partner mechanism [5l| for 
solving the doublet-triplet splitting problem, it is pos- 
sible to build models in which the proton decay rate is 
sufficiently suppressed [52J. An extended Higgs sector 
is required in that case. In principle, our model could 
be embedded in such extended SU(5) framework. For 
a review about proton stability, including also a section 
about SUSY SU(5) models, where further possibilities 
for sup pressing the proton decay are discussed, we refer 

to (HI . 



For what concerns leptogenesis, we argue that it can 
be realized in this model through the decay of the triplet 
or the singlet or both. However, since their exact masses 
are not determined from phenomenological constraints 
(contrary to the non-SUSY case addressed in Ref. 54]), 
it is not clear who is the responsible for leptogenesis: 
actually a combined action of the two could be possible, 
since their masses are of the same order of magnitude. 
To derive a definite conclusion as well as bounds on the 
parameters, a dedicated study would then be needed, 
which is beyond the scope of this work. 
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Appendix A: Renormalization group equations 

We present here the complete RGBs for this model. 
In the case of the MSSM parameters, we only explicitly 
write the new contributions to the 1-loop /3-functions, 
according to the definition: 



dt 



oMSSM.l o. 

Px + Px 



24,1 



(Al) 



where X can represent either Yukawa couplings, A-terms 



or soft masses. Clearly, all P^ ' vanishes below the scale 
of the involved 24 fields. The corresponding (3j^ ' can 
be found, for instance, in [53. For the new non-MSSM 
parameters, we provide the complete 1-loop /3-functions, 
still denoted as 13^' . 



J 



We first write the 



'X 



for the Yukawa couplings: 






(A 



VS J i 
24, 1\ 

VT )i 



^gl - igl + 3Tr(y„yJ) + 4 \ys\^ + ^ I^tI' + 3 \vv?)vs^ + \{vt- Vs) Vt. + iY,^Y:ys)^ 



-,g2-7.92 + 3Tr(y„yJ) 



l2/s|^ +4|yT|^ + 3\yvf ]yTi 



(2/I • yr) ys^ + {Y'eY:yT\ 



[PTv'), = ( - ^5? - igl - f 53' + 3Tr(y„y„t) + \y-'s\' + \ Wt? + 6 \yv\^)yv^ + 2iY,YJyy), 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
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The 1-loop /3-functions for the soft scalar masses read 
_ 3 

"I, ' ij 



{yk ivl K)o + (™i y*s)^ vl, + 2 yS. yl; ("^L + m|) + 2 ^j, a^^.) , (as) 

fc:') - 2(ye,, (yf m|J, + (m^-^ y^), y^^ + 2 y^, y^^- [m^^ + m|) + 2 A*y, A^^ , (A9) 

^™t' " 2(-4TO|yT + y^s^lvs + ■iylm^yv + "^H„(2l^^l' + '^^1' + 3|yy|2) + 

^m2,|y-r|2 + ™||rsP + 3m^|2;v|2 + ^\At? + jA^p + 3|/yp) , (AlO) 

P^ - 4(y]m|y5 + {ml^ + m|)|y-J|' + \As?) , (All) 

/3,^J^ - 2{ylmlyT + {m%^ + m|)|yT|' + \At\^) - IGM^g^ , (A12) 

P^J^ - 2(y^m|yv + (ml^ + ml.)\yv\' + \Av\') - yM^g^ _ QMlgl - ^M^gl - gjS , (A13) 

while {P^2 )ij = {P^2 )ij = (/3„2 )ij = 13^2 = 0. The hypercharge D-term contribution 5* is given by: 

5* = rriff^ — ra^Hd + Tr(— m| + m?e — 2m?c + m|^ + m|j) + 5(my — m|,) . (A14) 

Let us finally write the /3-functions for the trilinear terms: 

(/Sa!'')^ . = (^lyTp + \ys? + i\yv?) {Auh + (3yT • At + 2ys ■ As + &yv ■ Ay) (r„)y , (A15) 

(/?^^/) ^^ = 2yv^{ylAd), + ^Av^iylYa), , (A16) 

(/3't'') ^^. = \{AeyT)^yTJ + (^eyD^^s, + 3(r,y;), A?^- + 2{Y,y*s),Al^ , (A17) 

(/3lt'')^ = - l{As^ - 2^'hys^)gl - 3(^5. - 2M2ys^)gl + (3Tr(i;i;t) + h^^^\2 ^ ^j^^p ^ 3|^^,|2)^^^+ 

(6Tr(A„i;t) + 3^-^ . At + ITs ■ As + 6yv ' Ay) ys^+ 

(r/y;^), + 2{AJY:ys)^ + ^{y^ ■ ys)An + 3{yf ■ As)yT^ , (A18) 

(/^It'O, = " h'^T' " 2MiyT.).g? - 7(^T« - 2M2y„)52 + (3Tr(y„yJ) + ^lyrP + |ysl' + 3|yy|2) At.+ 

(6Tr(A„i;t) + ^y-^ . At + 2y* ■ As + 6y^ ■ Av)yn + 

{y:^y:At)^ + 2{AlY:yT\ + [y-^ ■ yT)As^ + 2{y* ■ AT)ysH , (A19) 

(/^It'O, = " T^(^^* " 2M^yv^)gl - 3{Av^ - 2M2yv^)gl - ^{Av^ - 2Nhyvi)gl+ 

(3Tr(i;,rJ) + ^lyTp + Ws? + ^\yv?)Av^ + (3Tr(A„i;t) + 3^-^ . At + 2y^ ■ As + lly^ ■ Av)yv^+ 

2iYdYJAv)^ + ^AdYJyv)^ . (A20) 
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